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We use a microscopic model, the Mirror-Oscillator-Field (MOF) model proposed in [T], to describe 
the quantum entanglement between a mirror’s center of mass (CoM) motion and a field. In contrast 
with the conventional approach where the mirror-field entanglement is understood as arising from 
the radiation pressure of an optical field inducing the motion of the mirror’s CoM, the MOF model 
incorporates the dynamics of the internal degrees of freedom of the mirror that couple to the optical 
field directly. The major advantage in this approach is that it provides a self-consistent treatment 
of the three pertinent subsystems (the mirror’s CoM motion, its internal degrees of freedom and 
the field) including their back-actions on each other, thereby giving a more accurate account of the 
quantum correlations between the individual subsystems. The optical and the mechanical properties 
of a mirror arising from its dynamical interaction with a quantum field are obtained without imposing 
any boundary conditions on the field additionally as is done in the conventional way. As one 
of the new physical features that arise from this self-consistent treatment of the coupled optics 
and mechanics behavior we observe a coherent transfer of quantum correlations from the field to 
the mirror via its internal degrees of freedom. We find the quantum entanglement between the 
optical field and the mirror’s center of mass motion upon coarse-graining over the internal degree 
of freedom. Further, we show that in certain parameter regimes the mirror-field entanglement is 
enhanced when the field interacts resonantly with the mirror’s internal degree of freedom, a new 
result which highlights the importance of including the internal structure of the mirror in quantum 
optomechanical considerations. 


I. INTRODUCTION 

Optomechanics describes the interaction of light with mechanical systems. When an optical field interacts with a 
mechanical object there is a redistribution of the photon momentum upon reflection. At the microscopic level this 
interaction results from the coupling of the electromagnetic (EM) field with the surface charges of the mechanical 
object (or the electrons in an atom). Still, the conventional approach towards studying optomechanical interactions 
only considers the effective boundary conditions for the optical field at the position of the mirror’s CoM that arise 
from the microscopic picture in the steady state limit. While the role of these internal degrees of freedom is universally 
acknowledged in the case of atom-field interactions when describing the mechanical effects of a field on an atom mm , 
their relevance in determining the optomechanical properties of larger systems is seldom discussed. 

However, questions pertaining to the transfer of quantum correlations, such as calculating the quantum entangle¬ 
ment between the mirror and the field [4-[10] or the superposition of two mirrors m, demand that we take into 
consideration the full quantum nature of the macroscopic object including the dynamics of its quantal internal degrees 
of freedom (similar to the two-level atom) to systematically account for all the quantum correlations present between 
the individual subsystems. Such a treatment becomes one of a practical necessity when studying the optomechanical 
entanglement for well-isolated systems that preserve coherences for longer time scales, for example when considering 
the quantum entanglement between the motion of atoms or atomic ensembles and a field. 

In a recent paper Galley, Behunin and one of the present authors [T] constructed a microphysics model called the 
mirror-oscillator-field (MOF) model that takes into consideration the microscopies of optomechanical interactions, 
providing a physically more complete theory for quantum optomechanics (QOM). The optical properties of the mirror 
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are captured in this model by introducing an internal degree of freedom of the mirror, referred to as the mirror- 
oscillator, or mirosc in [Tj. The mirror as an optomechanical element is described by two separate degrees of freedom 
corresponding to its center of mass motion (mechanics) and the surface charge that couples with the field (optics). We 
henceforth refer to these two degrees of freedom as the mechanical degree of freedom ( mdf) and the internal degree 
of freedom ( idf ) respectively. The idf and the mdf are each depicted by a quantum oscillator, with the idf coupled 
to an optical field that is modeled in (Tj by a massless scalar field. The idf is what provides the indirect interaction 
between the field and the mirror’s CoM motion, with its amplitude taking on field values at the position of the CoM. 
Compared to the traditional approach where the effect of the mirror on the field is represented by imposing boundary 
conditions on the field at the position of the mirror, a microscopic treatment captures the mirror-field interaction 
in a more physically consistent way as both the internal and mechanical degrees of freedom of the mirror enter in 
determining the dynamics self-consistently. As shown in [T], different parameters of the idf can describe a range of 
optical activities, from broadband to narrow band reflectivity. With specific parameter choices the authors in [T] 
made connections to well-known optomechanical models including those of Barton & Calogeracos El, Law m and 
Golestanian & Kardar |14j . 

The advantages of the MOF model over the usual practice of imposing boundary conditions and the role of the 
idf in capturing additional physical phenomenon is further expounded in this paper. We also use the MOF model 
to study the quantum entanglement between the mechanical motion of the mirror and the field. The conventional 
mechanism of mirror-field entanglement is by means of the radiation pressure coupling due to the photons of an optical 
field impinging on a mirror, transferring momentum to its center of mass [4]. There is virtually no consideration of 
how the mirror’s internal dynamics that give rise to its optical properties affect the entanglement of its external or 
mechanical degree of freedom with the field. Even though the full description of this interaction at the microscopic 
level is quite complex, to gain a qualitative understanding of the coupled interplay of the optical and mechanical 
degrees of freedom the present relatively simple MOF model can serve the purpose aptly and economically. As we 
shall see in this work, in some parameter regimes the dynamics of a mirror’s idf play a nontrivial and even a decisive 
role in determining the transfer of correlations and hence the entanglement between the mdf and the quantum field. 
We expect likewise for other related QOM effects. 

In the rest of this introduction we give a summary of the MOF model, followed by a description of the classical 
mechanical and optical properties of the MOF model in Sec. II. Sec. Ill treats the quantum dynamics of the three 
interlinked subsystems - the idf, the mdf and the field - which leads to all the interesting physical phenomena in 
QOM. In particular we show that the usual radiation pressure coupling is recovered as an approximation of the MOF 
model but one can go beyond these approximations to see new physical effects. With the solutions of the dynamical 
equations for this system we proceed in Sec. IV to derive the covariance matrix and calculate the entanglement 
between the mirror’s center of mass and the quantum field. The role played by the internal degree of freedom of the 
mirror is highlighted. We conclude in Sec. V with a discussion of the main points. 


The Mirror-Oscillator-Field (MOF) Model 


An optomechanical system consists of at least two components: a mirror interacting with a quantum field - where 
the "mechanics" refers to the mirror motion and "opto" refers to the field. One can think of three levels of description 
for this interaction: classical, semiclassical and quantum (see, e.g., m) in analogy to the studies of atom-field 
interaction. In the simplest classical electromagnetic description, the mirror is coupled to the field via the radiation 
pressure which can be obtained from the momentum flux imparted by the EM field that goes as the time-averaged 
Poynting vector of the field divided by the speed of light. In the semiclassical picture this force arises from the 
momentum transfer by the photons hitting the mirror, equivalently understood in terms of the gradient of energy 
density of the EM field on displacing the mirror. This leads to an interaction Hamiltonian of the intensity-position 
coupling form that goes as ~ Nx, where N denotes the photon number operator and x the displacement of the mirror. 
In a microscopic description of the radiation pressure force, one can think of the field discontinuity at the mirror’s 
CoM position inducing surface charge currents, which in turn experience a Lorentz force in the presence of the field 
that amounts to the radiation pressure force on the mirror’s CoM El- When one wishes to probe into issues like 
entanglement between a mirror and a field, a similar detailed treatment of the mirror and its interaction with the 
field is necessary, one that accounts for the transfer of quantum correlations between the mirror motion and the field 
as accurately as possible. 

Consider, for example, a single atom as an optomechanical element whose idf is represented by a two level system. 
We know that the interaction between the field and the atom’s two-level internal degree of freedom via photon 
emission and absorption is much stronger than the effective interaction of the field with the atom’s center of mass 
degree of freedom. The coupling between the optical field and atomic motion arises as the CoM motion alters the 
field configuration, thereby affecting the atom’s internal level activities. Thus when dealing with the case of atoms 
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Figure 1: Schematic representation of the interaction of a mirror with a field via its internal degree of freedom 


as an optomechanical element m-m , one needs to regard its internal level dynamics with careful consideration. 
Similarly, it has also been shown that the atom’s motional degree of freedom can affect the activities of its internal 
degrees of freedom such as spontaneous emission or motional decoherence as in f251E3] . For this reason, we point out 
the inadequacy in the study of mirror-field entanglement as one needs to take into consideration the idf of the mirror 
that is essential in mediating the quantum correlations between the optical field and the mechanical motion. This 
was one of the primary motivations in the construction of the MOF model in [I], which is highlighted in this work. 

Let us consider a point mirror interacting with a massless scalar field in (l+l)-dimensional space-time, the mirror 
is described by the two independent degrees of freedom - the mdf that has a mass M and is suspended in a harmonic 
potential of frequency 13 in addition to the idf described by another harmonic oscillator of mass m and frequency Q, 
as shown in Figffl 

While the mdf does not interact with the field itself, the idf is bilinearly coupled to the quantum field and 
constrained to be at the center of mass position leading to an effective interaction between the field and the mdf what 
we observe as the radiation pressure. The id/-field interaction determines all the optical properties of the mirror as 
has also been studied in pQ. We assume that the id/-field dynamics that represent the electronic excitations for the 
case of an atom happen at much faster time scales compared to those of the mechanical motion of the atomic center 
of mass, such that f l 13. 

For a non-relativistically moving mirror in the MOF model, the action is given by 

5 = J dt (QlkfZ 2 - ^M13 2 Z 2 ^J + (^mq 2 ~ \m^ 2 q 2 \ + j dz|((d t $ 2 ) - c 2 (d x $) 2 + \q$6(x - Z)^j (1) 


where we denote the center of mass position of the mdf by Z(t), the amplitude of the idf by q(t) and the scalar field 
by For relativistic motion which is required for the treatment of acceleration radiation such as the Unruh 

effect, one needs to use the proper time, modify the kinetic terms, and take care of the time-slicing scheme. Then 
the model will become a generalization of the Unruh-DeWitt detector theory |25l [25| . In drawing a correspondence 
between the scalar field and an electromagnetic field, we observe that the free field Lagrangian would correspond to 
that of an EM field if we choose Q(x,t) to represent the vector potential A. We have chosen a form of the bilinear 
coupling motivated by the electrodynamic form of interaction (~ — p ■ A), bearing in mind that the mirror’s idf can 
potentially represent the electronic level structure inside an atom. We note that this is different from the form of 
coupling in the original MOF model ID (~ Vi>). The 8(x — Z) factor in the coupling restricts the id/-field interaction 
to the center of mass position and the position dependence of the scalar field in turn leads to an effective force on the 
mdf We choose the coupling A to have the dimensions of the electronic charge e and t) to have the dimensions of 
A/c. This is in agreement with the correspondence of the MOF model with the Barton-Calogeracos (BC) model |12| . 
where in the limit of adiabatic idf evolution the coupling A can be physically identified as the surface charge density. 
Noticing that the free space permittivity in (l+l)-dimensions scales as eo ~ (Charge) 2 (Time)”(Mass) - (Length) - , 
the free field Lagrangian in jT] has been rescaled here by a factor of 6 q for dimensional consistency. 


II. CLASSICAL OPTOMECHANICAL PROPERTIES 

In this section we will illustrate how the MOF model can describe the classical optical and mechanical properties 
exhibited by a mirror, leading to the known intensity-position radiation pressure coupling. We begin with deriving 
the coupled equations of motion for the classical amplitudes of the mdf idf and field ({ Z , Z , q , q , <f>, $} respectively) 
from the action in ([l]) ( SS = 0) 
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Z + 15 2 Z= ^d x *(Z,t) 

(2) 

ii + n 2 q = --$(z,t) 

m 

(3) 

e 0 (x, t) — c 2 d 2 <h ( x , t) \ = A qS(x - Z) 

(4) 


It can be seen that the moving idf acts as a point source for the field and the idf is in turn driven by 4> at 
the center of mass position Z, which in the electromagnetic correspondence represents the electric field at the CoM 
position ~ E^j. Also, with A representing the charge density, it can be seen from ^ that the force on the surface 

charge degree of freedom goes as ~ A$. We have assumed here that the mirror center of mass velocity is in the 


non-relativistic limit, such that 


d z 

d t 


< c. For a relativistically moving mirror, the idf would more generally observe 


a Doppler shift of the field with respect to the moving center of mass as © becomes 


q + XI 2 q =- 

me 


Z°d, 


z 2 d x ) $(z") 


(5) 


where Z M Z 1 (t)) is the worldline of the mirror parametrized by its proper time r, and O = dO/dr. As the 

motion of the mirror center of mass leads to the motion of the charges sitting on the surface that interact with the 
field, the surface charges experience a Doppler shifted field which in turn changes their optical response leading to 
dynamically changing boundary conditions observed by the field. 

If one prefers to think in terms of applying boundary conditions on the field, in the MOF model it would correspond 
to the steady state response of the internal degree of freedom. Thus our model captures the full dynamical interplay 
as opposed to the static condition in the conventional approach. In fact, a simple generalization of the set up here can 
deal with a relativistically moving mirror as in dynamical Casimir effect (DCE), whereas the conventional method of 
imposing a static boundary condition on the field would fail to address dynamical situations wherein the time scales 
of the mechanical motion are comparable those of idf -field interaction dynamics. 

For now, we restrict our attention to a non-relativistically moving mirror. For the case where the system dynamics 
is driven by an incident field (and not by any external agent which accelerates the mirror as is in the setup of the 
Unruh effect), this is ensured from the separation of the timescales for the internal and center of mass degrees of 
freedom (f2 >> 13). We will demonstrate this further in Section IIB for the case of a single mode field. 

Knowing the coupled system dynamics, below we first look at how the radiation pressure force arises from our 
model in the non-relativistic limit. 


A. Classical Radiation Pressure Force 

As seen from ©, the mdf is driven non-linearly by both the idf and the field. We now eliminate the idf from the 
picture to obtain the mechanical force on the center of mass. 

From integrating the field equation of motion Q around the mirror center of mass position Z, we see that there is 
a discontinuity in the field spatial derivative. This can be understood as the discontinuity in the magnetic field across 
the mirror surface in the electromagnetic correspondence (d x <& ~ B) coming from the surface charge current q. In 
the non-relativistic limit we find the surface charge current as 

Xq = -e 0 c 2 5 a ,$||t (6) 

The surface charge current being induced by the discontinuous magnetic field across the center of mass position can be 
interpreted as the Ampere’s law in 1+1 dimensions. We eliminate the idf from the center of mass dynamics, defining 
the spatial derivative of the field at the center of mass position as d x 4> ( Z,t ) = (9 a -4> (Z + ,t) + d x <& ( Z~,t )) /2. We 
rewrite CoM dynamics as 

Z + V 2 Z = -- i - r e 0 c 2 (d x $) 2 \l + _ (7) 


1 This is an important point long explored and resolved in cosmological particle creation which results from the same mechanism but 
with the expanding universe playing the role of an external agent as in DCE. 
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We notice that the right hand side corresponds to the well-known radiation pressure force seen by a mirror in 

the non-relativistic limit m ■ This is justified based on the fact that the electric field vanishes at the mirror position 

n 2 

in the co-moving reference frame and the force being proportional to the EM field energy density then goes as ~ ^. 
To compare with the expression in m, we notice that for a perfect mirror there is no field energy density on one side 
of the mirror (d x <I> ( Z + ,t ) = 0) and we reduce to the known result. For an imperfect mirror there is a finite energy 
density of the EM field on either side of the surface, hence the net radiation pressure force is given by the difference 
in the field energy density on either side of the surface as in Q . In the MOF model the radiation pressure force can 
be interpreted as the Lorentz force arising from the interaction of the induced surface charge current ([ 6 ]) with the 
magnetic field Such an interpretation of the radiation pressure force as the Lorentz force on induced surface 

charge currents has been discussed in detail in M- 

Thus we have arrived at the classical radiation pressure force on the mirror in the non-relativistic CoM motion limit 
as one would find from imposing fixed boundary conditions on the field. Rather, in this case the boundary conditions 
resulting from the mirror-field coupling arise self consistently from the dynamical interaction between the moving idf 
and the field, as does the radiation pressure. We note here that while the form of the radiation pressure force we 
obtain from including the idf is identical to what we get from imposing the fixed boundary conditions, the boundary 
conditions themselves rather than being fixed are determined by the dynamics of the id/-field interaction. This more 
generally includes retarded influence of the moving surface charges on the field in a dynamical way. To see this more 
concretely, consider the idf amplitude solution from ©> 

q{t) = Qh(t) + I dt'Gift -1') 

Jo 


- d>{Z(t'),t') 


( 8 ) 


where we define qh as the homogeneous solution for the free idf evolution and Gift — t') = * )"> as ^h e Green’s 

function for the idf. We use this to eliminate the idf from the field’s equation of motion to get 


\2 _ ft 

e 0 (<9 2 - c 2 d 2 ) $ (x, t ) -I- S (x - Z(t]) / d t'd t Gi (t - t') ${Z{t'), t') = Xqh.S (x - Z(t )) (9) 

m Jo 

We see that the idf is driven by the field and influences the field in return, as captured in the second term on the 
left hand side that represents the retarded influence of the idf on the field, meaning that the radiation pressure force 
depends on the coupled non-Markovian dynamics of the field, center of mass and the idf. Thus, we can identify the 
term dtGi (t — t') = x (t — t') as the susceptibility function for the mirror. To compare with the case where one applies 
boundary conditions as opposed to including the idf dynamics self-consistently one needs to include the coupling of 
the idf with a bath so as to reach the steady state response of the damped idf. We will further illustrate this point 
and the role of the internal degree of freedom in determining the optical properties of the mirror in the following 
subsection. 


B. Optical properties 

To study the optical properties arising from the MOF model let us consider a single mode field at frequency w and 
amplitude $0 driving the mirror’s idf. Assuming that the mirror’s CoM is at the origin in equilibrium, we make the 
following plane-wave ansatz for the field 

f x t ) = (&(—x) (e ikx + R(Lo)e~ ikx ) + 0(x)T(u})e ikx ) + H.C. (10) 

U) 

where we have introduced the frequency normalization factor (Sl/ui) to take care of the fact that in the EM corre¬ 
spondence the electric field amplitude (E ~ d t A) is independent of the frequency of the field. R (w) and T (oj) refer 
to the reflection and transmission coefficients of the point mirror, such that T (oj) = 1 + R(oj). In considering the 
interaction of the idf with only a single field mode, we include a damping (7 f) that arises from its coupling with the 
remaining field modes. For current purposes, we assume that the damping is small ( 7 f <C f2) so that one can ignore 
the dissipation of the incident plane wave. As in [T|, we assume that in the steady state regime the idf oscillates at 
the frequency of the incident field. In which case, we find 

q(t) = r T A o + H.C. 

m (oj z — U z ) uj 


(ii) 
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From the mirror center of mass dynamics ([t]) , we can see that in the presence of the incident drive the center of 
mass consists of a time-independent and a high frequency (2 ui) radiation pressure term, coming from the non-linear 
interaction of the incident B field and the induced surface charge current. In the limit 13 <C f2, the high frequency 

component of the mirror amplitude denoted by Z% u scales as jZ^j ~ > which, in the near field- idf resonance 

regime (to ~ f2), is much smaller compared with the mirror amplitude coming from the constant radiation pressure 

part Z 0 ~ e °Miy2 , noting that ~ 02 1- Thus we find that the mirror position evolves essentially at its natural 

frequency 13 under the constant force. 

We assume that at the classical level the center of mass motion does not affect the id/-held coupling and the 
resulting optical properties from the interaction. More explicitly, the phase of the held mode that is resonant with the 
idf changes by a very small amount over the length scales of one amplitude of the mdf that is A if = (f l/c) Z 0 <C 1. 
This restricts the held amplitude to 


_ |2 M13 2 c 

'•’■"I « w 


( 12 ) 


This is a self-consistent validity constraint which ensures that the optical properties of the mirror are unaffected by 
the center of mass motion to hrst order, to reaffirm our plane wave ansatz (101. Physically speaking we assert that 


the mirror CoM motion is much smaller than the wavelengths of the held that it interacts with. The sub-wavelength 
motion approximation is valid for the case of trapped atoms spatially confined in a harmonic trap (trap frequency 
being 13 in this case), interacting with an optical held of frequency u>. 

In the plane wave ansatz, we find the surface charge current for the idf © as 


A q = —eoc 2 d x & (x, t) Ift sa — 2ifceoc 2 $ 0 —e lult R (w) + H.C. 

CO 


a -iuit ] 
—iujt ; 


= -2 ie 0 nc<3> 0 e- lult R (w) + H.C. 


(13) 


We can notice here that the induced surface charge current is proportional to the mirror rehectivity. Thus, as expected, 
a higher rehectivity leads to a larger radiation pressure force. 

Now within the non-relativistic and sub-wavelength CoM motion approximations, we consider the MOF model with 
the two different forms for the coupling term - (1) < 7 <I> (as previously analyzed in |T|) and (2) - and study the optical 

properties that arise from these two couplings in different parameter regimes. 


1. g<f> coupling 


Let us hrst consider the coupling as in pQ and start with drawing the correspondence between the interaction 
term for the scalar held vis-a-vis an EM held. As motivated in the section II.B.l in ||I| when comparing the MOF 
model with the Barton-Calogeracos (BC) model, we choose the coupling A to have dimensions of the charge density 
such that dimensionally A ~ (Charge) (Length) -1 . Going back to the interaction term in the original action we use 
this to find the dimensions of the scalar held as <1> ~ (Mass)(Length) (Time) - /(Charge) and rescale the free held 
term accordingly, we get for the free held action 


S F = 


ep 

2c 2 


d t /dx((d t $) 2 -c 2 (d^) 2 ) 


This leads to the coupled idf held equations of motion for a hxed center of mass as 

e 0 /c 2 (5 t 2 $ - c 2 d 2 <h) = Xq6(x) 
mq + mSfq = A<h(0, t) 


(14) 


(15) 

(16) 


For a plane wave incident on the mirror, using the ansatz (101 to solve for the surface charge current as in (13) in 


the steady state limit we get the rehectivity for the case of q& coupling as 

—iA 2 c 




iA 2 c + 2mujeo(uj 2 — f l 2 )' 


\R\ 2 = 


1 + r 2 ?7 2 (1 — if) 


(17) 


where we have defined the ratio of the held to the idf frequency as r] = ui/Cl and r p = ^ 2 ™? e e ° ^ = Q,/flp. We identify 
\ 2 

the quantity f Ip = 2 mU^e as plasma frequency, again motivated by the comparison with the BC model. As found 
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in jT], the mirror becomes perfectly reflecting for 1) infinitely strong idf -field coupling, A — > oo, 2) perfect resonance 
between the idf and the incident field, w = fl or 3) massless idf m —> 0. Now observing that the reflection spectrum 
is completely characterized by the two frequency ratios r p (ratio of the idf to plasma frequency) and r] (ratio of the 
field to idf frequency), we consider different values for the parameter r p and look at the reflectance as a function of 
the field frequency for a fixed plasma frequency, as shown in Fig. [2j 

To invoke the correspondence with the BC model m we need to assume that the idf evolves adiabatically in the 
limit {to —>• 0, ST —>• oo} such that the quantity mCl 2 = n that physically corresponds to the mass density of the surface 
charges stays finite. In this limit since r p 3> 1 (fl —t oo), we see a resonant behavior in the reflection spectrum around 
the idf frequency Cl. In the regime where r p <C 1, the reflection spectrum shows a high frequency cutoff behavior 
similar to the case of bulk metals with Drude-Lorentz response. As shown in Fig[2j given the plasma frequency for 
silver (Clp = 1.37 x 10 16 Hz), we compare the known optical response with our model and find that a idf to plasma 
frequency ratio r p ~ 0.3 mimics the cut-off behavior reasonably well. Knowing that the charge carrier density for 
silver is n s = 5.8 x 10 28 to - 3 and using the BC correspondence to find A = n s e, we can deduce all three idf parameter 
values. 




1.0 £. £..0 O 0.£. O.H 0.0 0.0 I ! .£. I .H / .O 1.0 O 

2ttC/uj xIO' 6 2 tt dto xIO' 7 


(a) 


(b) 


Figure 2: Reflection properties from the two different forms of coupling (g$ and g<l>) (a) Reflectance as a function 
the incident field wavelength for different idf to plasma frequency ratios ( r p = Cl/Clp), the plasma frequency is fixed 
at Clp = 1.37 x 10 16 Hz (for silver) from gd> coupling, choosing r p ss 0.3 mimics the cut-off behavior for silver (b) 
Reflectance and transmittance spectrum from g<I> coupling to simulate the optical response for a photonic crystal as 
from the experimental results in m Each resonance corresponds to a separate idf. 


2. g<f> coupling 

As in the previous subsection we find the reflection coefficient for the g>I> coupling as 

—i\ 2 u> 2 to 2 /Cl 2 r/ 2 

R(u) = -X--yyyy; |R(w)| = --yfy-yry (18) 

iA-w + 2me 0 c(uj 2 - Cl 2 ) ui 2 /Cl 2 + (oj 2 /Cl 2 - l) 2 + r 2 (?? 2 - l) 2 

where again we have defined the ratio of the field to idf frequency as rj = ui/Cl and the ratio of idf to plasma frequency 
as r p = Cl/Clp = redefining the plasma frequency as Clp = 2T ^ eoC - We see that mirror becomes perfectly 

reflecting for the same conditions as in the case of coupling {A —>• oo,co = Cl,m —> 0}. Unlike the gT coupling, we 
do not see a perfect reflection at ui = 0 which was an artifact of the monopole coupling between the idf and the field. 

The optical response exhibits a resonant behavior around the idf frequency Cl, since the reflectivity is maximum 
for Tj = 1. For this reason it is natural to consider optomechanical elements with built in resonances such as photonic 
crystals or atoms as an application. We find that one can mimic the optical response of a photonic crystal structure 
(see Fig{2| by choosing the resonant frequency of the idf as the resonant mode of the photonic crystal, for multiple 
resonances we choose multiple internal degrees of freedom such that Cli = w} es at each resonance peak. The sharpness 
of the resonance is determined by the quantity r p , since the parameter r p determines the coupling strength of the field 
to a particular resonance mode of the structure. Thus one can determine the two parameters that characterize the 



































optical response in our model, namely r p and f 1. To completely determine all the parameters of the internal degree 
of freedom {m, Cl, A} we need to draw an additional physical correspondence between the internal degree of freedom 
and the physical setup as we did for the previous case of g<E> coupling by identifying the coupling constant A as the 
charge density. 

As we had noticed previously, the mirror reflectivity characterized by the idf parameters determines the strength 
of the induced surface charge current (13I which in turn factors into determining the radiation pressure coupling. In 
the following section we will show that the same applies to the case of coupling between the quantum fluctuations of 
the mirror and the held. We now turn to look at the coupled quantum dynamics of the three subsystems in the MOF 
model. 


III. QUANTUM DYNAMICS OF THE COUPLED MIRROR-OSCILLATOR-FIELD (MOF) SYSTEM 


Let us perturb the original action about the classical solutions as {Z + Z,q + < 7 , $ + $}, with O being the 
deviations about the classical solutions O. Assuming that the center of mass motion about Z is small and restricted 
to the sub-wavelength regimes (kZ <C 1 ^ for the held modes below a certain high frequency cutoff, we expand the 
action up to third order in the fluctuations about the classical solutions ignoring terms that are second order or higher 
in kZ. We go up to third order to specifically include the term that couples the perturbations of all three subsystems 
(labeled as MOF below) to arrive at the non-linear intensity-position Nx'j coupling. In the subsequent dynamics 

we shall only consider bilinear interaction terms to preserve Gaussianity of the individual subsystems. As we will see, 
truncating the action up to second order corresponds to the linearized approximation in the limit of strong mean-held 
amplitude, also called a background held expansion in held theory. 

We write the perturbed action as 


S 3 = / dt ( i MZ 2 - ^MU 2 Z 2 J + Qrrfy 2 - ^mCl 2 q 2 


OF 


mdf (M) 


da; 


idf (O) 


(a t $) 2 -c 2 (s^) 2 ] +A5 (x-Z) 


Field (F) 


q$ +q (d x $) Z + f {d x §) Z + f (d„<I>) Z 


MF 


OM 


MOF 


(19) 


One can observe several points from the above expression, hrstly, we find that there is an effective coupling between 
the fluctuations of the mirror center of mass and the held via the internal degree of freedom as denoted by the terms 
MF and MOF. To the lowest order, the mirror-held coupling strength is proportional to the classical surface current 
q , implying that the fluctuations of the held are the most sensitive to the huctuations of the mirror center of mass if 
the surface current is at its largest. In the single held mode case this is proportional to the reflection coefficient of the 
mirror as seen in (13), meaning that a highly reflecting mirror leads to large effective MF coupling strength. Secondly, 


there is also an effective coupling between the idf and the mdf huctuations denoted by the terms OM and MOF, 
which to the lowest order is proportional to the spatial derivative of the held (or B) at the center of mass position. 
The coupling strengths of the interaction terms between the idf and the mirror (OM), and the held and the mirror 
(MF) are determined by the classical solutions of the held and idf amplitudes as found in the previous sections. 

We get the following equations of motion for the coupled mirror and held dynamics 


Z + 13 2 Z = 


A 
M L 


qdji> (Z, t) + q (Z, t) + d x t> [Z, t ) j 


e 0 (d 2 <& - c 2 d 2 ^j = XqS (x - Z) - A (q + q) d x ( 6 ( x - Z)) Z, 



( 20 ) 

>-z))z 

( 21 ) 


It can be seen here that unlike the classical equations of motion, the held huctuations are not only driven by the idf 
but also by the huctuations of the center of mass position. From integrating (21) around the classical center of mass 
position Z , we get the surface current fluctuation as 


A q = -e 0 c 2 d x <&\ 


(22) 
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just as the classical version interpreted as the Ampere’s Law in 1+1 D in Q. Using this and the classical surface 
current to eliminate the idf from the center of mass dynamics (20), we get 


Z + 13 2 Z = 


-eo c 
M 


(a x $) (d x ®)\l + + l (M) 


(23) 


We can see that the first term on the right side corresponds to the radiation pressure coupling in the linearized 
approximation which is valid for large photon numbers in the presence of a classically driven field. The second 
term goes beyond this approximation, which corresponds to the Nx type of coupling, required for treating situations 
with small photon numbers. Considering $ represents the quantum fluctuations of the field, we can understand the 
radiation pressure force at the quantum level as arising from the asymmetry in the field fluctuations on either side 
of the mirror. Say, if there were a cavity present on one side and free space on the other, the radiation force from 
the cavity side would be stronger in comparison because of the small quantization volume leading to asymmetry in 
the density of field modes as in the case of Casimir force [28] ■ Such an interpretation of Casimir force as a radiation 
pressure force from the vacuum field has been discussed by Milonni et al in [29] for the case of two perfectly conducting 
plates. 

We now restrict ourselves to second order perturbations in the original action, to keep all the interaction terms 
bilinear such that starting out with Gaussian initial states for the three subsystems, Gaussianity of the individual 
subsystems is preserved. We derive the conjugate momenta from the second order action as 

p = mf + \$(Z,t) + \d x $(Z,t)Z (24) 

P = MZ (25) 

= e 0 $(x,t) (26) 


It can be seen that the fluctuations in the idf are influenced by both the mdf and the field and hence mediate the 
effective interactions between the two. Identifying the dynamical variables {Z,q, $} as the quantum fluctuations of 
the mdf, idf and the field respectively about their mean-field amplitudes, we arrive at the second order Hamiltonian 


= 


2 M 


+ 1 1,2 + P 2 + + + 


■)2 ~2 


2m 


mdf(M) 


d * ( 5T 0 + W 


(cU) 


idf(O) 

2 1 + ^<1 (z, t ) 2 - -mz, t ) - -a MZ, t)pz + —d x $(z, twz, t)z - xfdMz, t)z 


2m 


m 


m 


m 


Field (F) 


OF 


OM 


MF 


(27) 


We notice that the mdf now observes a renormalized oscillation frequency and the scalar field sees a shift coming from 
the term quadratic in $ which is analogous to the diamagnetic term ~ 2 e 2 A Z of the minimal coupling Hamiltonian. 
The bilinear interaction terms represent the coupling between the idf and the field (OF), mirror and the idf (OM) 
and mirror and the field (MF) respectively. Physically, the terms that are second order in A arise from the field-field, 
mirror-mirror and field-mirror couplings mediated via the quantum fluctuations of the idf. The terms that are first 
order in A in the couplings between the idf-mdf ( OM) and mdf -field (MF) fluctuations come from the classically driven 
solutions for the field and the idf respectively. Specifically, we note that the MF interaction contains two terms, the 
first one of which represents the effective mirror-field interaction mediated via the quantum fluctuations of the idf 
while the second one represents that from the classical surface charge currents. Since the conventional approach does 
not include the fluctuations of this extra quantum degree of freedom, it misses out on the fluctuation-mediated part 
of the effective mirror-field coupling. As we shall see later, this term becomes dominant in the strong coupling regime. 

We also note that in the absence of a classical drive, the only interaction is between the idf and the field (OF) up 
to second order. To be able to see an effective mirror-field interaction one needs to include third order terms in the 
fluctuations as illustrated before. 

In the following section we study the above Hamiltonian for the case of a driven single field mode and find the 
subsequent entanglement dynamics for the mirror CoM and the field, coarse-graining over the internal degrees of 
freedom. 
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IV. MIRROR-FIELD ENTANGLEMENT IN THE MOF MODEL 


Entanglement between a field and a mechanical oscillator has been widely studied in cavity optomechanical setups 
in several contexts m-m, with the essential mirror-field coupling mechanism being the radiation pressure wherein 
the field exerts a force on the mirror center of mass by means of photon-momentum transfer and observes a phase shift 
proportional to the mirror displacement in turn. We now look at the entanglement generation from a microscopic 
perspective as described by the MOF model, considering only a single mode of the scalar field in our model as in the 
usual cavity optomechanical setups to deduce some key physical features of the mirror-field entanglement that arise 
from the inclusion of the idf. 


We first simplify the Hamiltonian (27) for the case of a single field mode that is being externally driven to look at the 


dynamics of the coupled MOF system and then coarse-grain the idf to find the sought after mirror-field entanglement. 
Consider the scalar field in a region of length L (assuming L approaches infinity), the field fluctuations can then be 

„ik r , 


2 injfi e o L 


[a n e 


ale~ ik - 


written as the sum of all discrete modes of the cavity of length L as $ (x,t) = 

with ajj and a n representing the creation and annihilation operators for the n th field mode. We pick a single field 
mode at frequency ui interacting with the point mirror at the origin (Z = 0) assuming that the center of mass motion 
is in the sub-wavelength regime as before. 


** (X, t ) = (a u e ikx + ~ale~ ikx ) 


(28) 


The above expression represents the fluctuations of the free field without any imposed boundary conditions unlike 
the standard treatment where the quantum fluctuations follow the mode functions of the classical field (see [361 for 
example). In the steady state, the strength of the field fluctuations would be determined by the boundary conditions 
as they emerge from the id/-field interaction self-consistently. 


For a single field mode, we rewrite the free Hamiltonian part in (27) as 


kfree + \ Ml5 ' 2 z 2 + m + 1) 


A 2 


2mojeoL 


Q'CjQ'U) 


A 2 


AmueoL 


((a u ) 2 + (at) 2 ) (29) 


Hm 


Ho 


Hf 


where we have redefined the dynamical variables associated with the idf in terms of the creation annihilation operators 
{, 6} as q = yJ (b + b^) and p = —iyj^^{b — b t^. The renormalized mechanical frequency is defined as 

O' 2 = 13 2 + (d x $(Z, t)) 2 . As mentioned in the previous section, the correction term (d x $(Z ,t))~) contains 

two contributions - a time dependent part oscillating at a frequency ~ 2 u> and a time-independent part. In the 
rotating wave approximation (RWA) the time dependent term can be neglected. However, if the field mode was 
resonant with the mdf one would see parametric amplification of the mirror center of mass motion m ■ For the free 
field part we notice that the interaction leads to an energy correction u> —> uj + A 2 / (2mojeoL) that is second order in 

A, this corresponds to the shift coming from the diamagnetic contribution 2 mc 2 ^ 2 ) f° r the EM case as indicated 

in the previous section. This also leads to the fast oscillating terms for the free field (~ 2w), which correspond to the 
photon-pair production and annihilation as in the case of dynamical Casimir effect |.'U)il.'L‘>j . Moving to the interaction 
picture with respect to Hq = Ho + Hp to eliminate the fast dynamics of the system and invoking RWA, we write the 
interaction Hamiltonian in a simplified form as 

H int = h (a OF tfae~ iAt + a OF ba f e iAt ) + H [a OM be iAt + a* OM b f e~ iAt ) (c + c f ) +H ( a MF a + a* MF a f ) (c + c f ) 

(30) 

Here we have defined the operators in the interaction picture as {a, a^} = {d u e lult , dte~ lujt } and {b, ih} = 
{&e* nt , b^e~ int } and the detuning A = w — H represents the detuning between the field and the idf. The oper¬ 
ators {c, c^} correspond to the creation and annihilation operators for the phononic excitations of the mdf, with 

Z = \J 2 mw + cl) = yj jjjyZ and P = — if bMpL (g _ gt) = y/fiMU'P. The operators Z and P are the dimen¬ 
sionless position and momentum fluctuations for the mirror center of mass. In moving to the interaction picture we 
have ignored the second order correction terms A 2 /to) in the free field Hamiltonian H F . 

The coefficients o^s represent the effective bilinear coupling strengths between the single excitations of the three 
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subsystems with 


ocof = — 


aoM = 


OiMF = 


iX / n 
2 V rnueoL 
S2<F 0 A 


Q 


mMl5' 


2c 

2c V MV'ue 0 L 


1 


iA 2 

-h 2eo cujR*(cj) 

m 


OAn 


h 


iX 2 


2MIX' V 2mcujeo 


+ R* (w) 


(31) 

(32) 

(33) 


where we have defined the dimensionless field amplitude Aq = <E>o/ \] 2 Jc. 0 L ' ^ can then be seen from (331 that for 

a perfectly reflecting mirror (R*(ui) —» — 1) the second term in the effective coupling strength cxmf between the mdf 
and the field is the same as what one finds from standard boundary condition approach (see Appendix |B| . 

We note that all the effective coupling strengths contain the idf mass and charge parameters in the combination 
~ A 2 /to which corresponds to the plasma frequency flp(= A 2 /(2tococ)), meaning that one can deduce all the effective 
single excitation couplings (a tJ s) from the two parameters that also completely characterize the reflection spectrum, f 1 
and flp as defined in section ['ll B| Thus given the reflection spectrum of a mirror, one can find the parameters f Ip and 
Q, knowing which one arrives at the various effective coupling strengths. Fig[3]shows the dependence of the reflection 
coefficient and these effective couplings on the dimensionless plasma frequency (Qp/Sl) and detuning (A/f2). 

It can also be observed from (31) (331 that the coupling strengths increase as the original idf- field coupling A 
increases and decrease as the idf mass to increases, meaning that a "lighter" idf leads to stronger effective coupling 
strengths. Also, a heavier mirror CoM couples more weakly to the idf and the field. The effective id/-field coupling 
aoF is independent of the driving field amplitude <f>o as expected, since as one sets the drive amplitude to zero it can 
be seen that there is no mirror-field and id/-mirror interaction in second order except the idf- field coupling from the 
direct interaction. To be able to see any coupling and hence entanglement between the mirror and the field in that 
case one needs to include the higher order terms as was discussed before. 

We also note here that in the weak coupling regime where flp <C 1, the mirror reflectivity and the effective mirror- 
field coupling strength olmf as a function of the idf-Held detuning peaks sharply at resonance (A = 0) as seen from 
Figj3ja) and Fig j3] (d). The field amplitude and detuning with respect to the idf change the coupling strengths 
appreciably. While in the standard treatment of mirror-field interactions via boundary conditions it is common to 
study the effect of the field intensity on the mirror-field coupling, we highlight that including the presence of idf lets 
us see the effect of the field-id/ detuning on the mirror-field interaction, allowing us to probe the effective coupling 
strength as a function of the reflection properties of the mirror. 

As noted before, the two terms in the effective mirror-field coupling olmf denote the interaction mediated via the 
quantum fluctuations of the idf and its classical amplitude respectively. The strong-coupling limit, where one would 
expect to see non-Markovian dynamics is also where the contribution from the idf fluctuations becomes substantial, 
as seen from the first term in (331. 


Now we use the interaction Hamiltonian (301 to write the equations of motion in terms of the coupling constants 
oto s as 


dZ 
dt 
d P 
d t 
dq 
di 
dp 
df 
d«F 
d t 

dn 

dt 


=13'P 

= — 15'Z - 2 (ReaoA/q — Imaojup) — 2 (R ectMF& - ImaMfll) — 7 P + £ 
=A p — — 21m olomZ 

= — A q — |a 0 F|n — 2ReaoM-^ 

=\o:of\q — 21 m c*mfZ 
=\&of\p — 2R eotMFZ 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


wherein we have redefined the slow moving dimensionless idf and the field quadratures as q = ——-, p = 

— i ——-, $ = a +g T and II = —i ■ Also, to account for the fluctuation-dissipation mechanism for the 

mirror center of mass resulting from its coupling to the thermal bath, we have introduced the mechanical damping 7 
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A/S! 


(a) 



A/a 


(b) 



A/a A/a 


(c) 


(d) 


Figure 3: (a) Reflectance as a function of the dimensionless parameter Up/U (ratio of the plasma frequency to the 
idf s natural frequency) and the idf- field detuning A/a. It can be seen that for weaker coupling corresponding to 
Up/U < 1, the reflection spectrum has a sharper resonance. The effective bilinear coupling strengths for both 
(b) idf-field («of) and (c )idf-mdf (qsom) increase with increasing plasma frequency as ~ \/Up. (d)The effective 
mdf- field coupling coefficient (cimf) hr the weak coupling limit is largely determined by the reflection coefficient, 
while for strong coupling the fluctuation mediated part becomes relevant, as can be seen from (33). 


and noise £ for the mirror. In the high temperature limit, the correlation function of the noise is given as 
2M^/kBT8{t — t'), with T as the temperature of the thermal bath. 


(K*)l(0) = 


At this point one can make a crucial observation that if we assume weak-coupling such that U p <C a, the fluctuation 

mediated part of the effective mirror-field coupling is negligible or ~ ( w )- Additionally, if we 

consider the field detuning to be small enough such that A <C U p , the reflection coefficient R*{co) ~ 1 leading to a 
nearly perfectly reflecting mirror. In this parameter regime, it can be seen that the mirror-field dynamics from (34)- 

f found from the MOF model reduces to that from the conventional boundary condition approach (see Appendix 
>r further details), provided that the idf contribution is negligibly small. To ensure which we consider that the 
idf is coupled to the continuum of field modes with a coupling of the form q&i, where represents the \ th field 
mode, leading to a damping coefficient 7 f. Also, to mimic the scattering of surface charges by lattice ions of the 
mirror, we introduce a dissipative bath of internal degrees of freedom such that each bath oscillator is coupled to the 
idf with a coupling of the form q ■ qi, where qi represents the position variable for the i th bath oscillator, giving an 
effective damping coefficient of 7 ^ for the idf. Using separation of time scales, we find the steady state idf amplitudes 


as q st = — and p st = , where the operators C/s stand for C\ = \cxof\& + ZhnaoM Z + £7 and 

C 2 = |«OF|n + 2ReaoM-^ + £U Now for the case of near perfect reflection since the detuning A is small, for the 
steady state amplitudes to vanish, we must have 7 ij A. It can then be seen that the dynamics obtained from the 
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two approaches agree perfectly with each other as shown in Fig. [4j 



Figure 4: Evolution of mirror-field entanglement as measured by the logarithmic negativity Emf (see Appendix [a] 
for definition) as obtained from the boundary condition approach (Appendix [B| and the coupled MOF dynamics. 

We find that for an isolated idf the time scale for entanglement is largely determined by the effective idf-Aeld 
coupling (aof). The two approaches concur in the weak coupling limit for a strongly damped idf. The parameters 
values, in units where c=l, h = 1, used here are m = 0.001, f l = 100, M = 10, 0 = 0.1, flp = 5, Aq = 10 -4 and 
T = 1000. The effective id/-field coupling strength |ckof|/G ~ 16. 


Now going back to the case of an isolated idf, we can identify the radiation pressure force from (351 as 

Frad = —2 (ReaoArQ — ImaoMP) — 2 (ReaMF^ — IiuomfII) 


(40) 


We can see that the linearized radiation pressure force depends on both the fluctuations of the idf and the field 
variables and is generally dependent on the idf parameters. Hence as long as the idf fluctuations are non-vanishing, 
the radiation pressure shot noise is determined by the shot noise of both the field and the idf, meaning that in order 
to go below the standard quantum limit for the radiation pressure force one needs to take into consideration the 
squeezing of the idf quadratures in addition to those of the field |341 [35]. 


From (34) (391 we write the solutions to the equations of motion for the idf and the field variables as 
q(t) = q h{t) + f d t'G 0 (t — t') (—A|aoF|n(f') — A aowiZ(t') + 2|a!0p|Ima! M pZ(f , )) 

Jo 

P (t) = p h (t) + [ d t'G 0 (t - t') (A| a 0 F\$(t') + 2|aoF|Rea M FZ(t') - aojuP(O) 

Jo 

*(i) = *„(t)+ [ dt , G / (t-t , )(A|aoF|p(f , )-2Ima M pP(t')) 

Jo 

II(t) = + f d t'Gf(t — t') (—Alaoplq)! 7 ) — \aoF\oiOM^{f) — 2ReaMFP(f , )) 

Jo 


(41) 

(42) 

(43) 

(44) 


Here we have defined the idf and the field Green’s functions as Gott) = sln( Vl ao - F h+ A anc j Gf(t) = sln (l a Q J | , l*) anc i 

J K ’ y|a 0 F| 2 +A 2 ’ I“Of| 

the homogeneous solutions as {q^, p / l7 11^}. It can be seen that the frequency of oscillations for the slow moving 
idf variables is Qidf = \/Wq F \ + A 2 and that for the slow moving field variables is H/ = |oof|- In the steady state 
limit, we can use these solutions to rewrite the equation of motion for the late time mirror CoM dynamics as 

M (<9 t 2 + O' 2 ) Z (t) + r y dZ ^ + d t’Goit - t') (-A a 2 OM + 2a 0 F^a M Fa 0 M) Z (t 1 ) = 

I &t'G 0 (t — t') (2aoM A|aoF|n(t')) + f At'Gf{t — t') (2ReaMF^\o:oF\p{t l ) + 2 ImaMFA|aoF|q(^ , )) + £ (45) 


On the left side one can identify the two terms in the integral as the retarded influence of mirror-idf-idf-mirror 
interaction and the mirror-idf-field-mirror interaction respectively. The first term on the right side denotes the mirror 
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being driven by the idf influenced field and the second term stands for the mirror being driven by the field-influenced 
idf. We can see that in the absence of any detuning the CoM motion is only driven by the thermal noise term. 

One can find the normal modes of the system from (34 ( (|39|) and their time evolution to obtain the 6 x 6 dimensional 


covariance matrix of the coupled MOF system numerically. We define the MOF covariance matrix Vmof as 


( Vmm Vmf Vom \ 
Vmf v ff V 0 f 
VSm v Sf Voo ) 


(46) 


where the on-diagonal sub-matrix Vkk stands for the covariance matrix of the k th reduced subsystem, defined as 
(Vkk)y = \ X^' 1 }^>, with xj k ^ and X^' 1 representing the i and j quadratures corresponding to the position 

and momentum variables of the k th reduced subsystem, more explicitly X^ p^}. Here, { 0 i, 02 } denotes 

the anti-commutator between the operators 0\ and Oi- The off-diagonal sub-matrix Vki consists of the correlations 
between the k th and the ^subsystems, such that (Vki), :j - = \ X^}^, where the i and j quadrature components 

belong to different subsystems. 

We then choose to look at the part of the covariance matrix that represents the mirror and field reduced covariance 
matrices and correlations, that is, 


Vmf 


V m 

V T 

v MF 


Vmf A 
V f J 


(47) 


and find the logarithmic negativity E^ F as obtained based on the positive partial transpose (PPT) criteria for 
determining separability (see Appendix [A] for details). It can be shown that calculating the MF entanglement from 
the sub covariance matrix Vmf is equivalent to coarse-graining over the internal degree of freedom and then finding 
the MF entanglement. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 


Sit x10 4 

Figure 5: Mirror-Field entanglement given by the logarithmic negativity for an undamped idf as a function of the 
dimensionless id/-field detuning (A/O) and dimensionless time (f H). We observe that the entanglement peaks for a 
resonant idf field interaction at ( A/13 = 0) and the oscillation time scales are determined by the effective idf- field 
coupling (|qiof|)- The parameters values, in units where c=l, H = 1, used here are m = 0.001, f 1 = 100, M = 10, 

0 = 0.1, n P = 0.05, A 0 = 10~ 4 and T = 1000. 

As was discussed before, at the id/-field resonance (A —» 0) the reflection coefficient and hence the effective mirror- 
field coupling strength goes to its maximum value (See Fig{3](a) and Fig. [3](d)) . As a result, we observe in Figj5]that 
there is a peak in the mirror-field entanglement near id/-field resonance. As was emphasized before, this effect is not 
considered in the standard treatment of optomechanical interactions since the internal degree is coarse-grained over a 
priori to arrive at the boundary conditions for the field. This effect is more pronounced in the weak-coupling regime 
where the reflection coefficient has a sharper peak at resonance as seen from Fig j3ja) and Fig(3jd) . It can also be 
observed that for A/O = —1 or equivalently Q = ui + 0, the entanglement is sustained for longer times. Physically, 
this pertains to the process wherein a field photon and a mirror phonon combine to give a single idf excitation (or vice 
versa), corresponding to the two-mode squeezing Hamiltonian which then entangles the field and the mirror modes 










15 


as a result of the interaction. Such an observation had also been made in []J| for the case of a cavity driven with a 
red detuned drive in the sideband resolved regime where it was shown that the steady state entanglement goes to 
a maximum when the cavity-drive detuning was equal to the mechanical oscillation frequency. Drawing an analogy 
between the two cases, we find that cavity resonance for the usual cavity optomechanical setups is similar to the idf 
in the MOF model in that they both mediate the interaction between the mechanical motion of the mirror and the 
external field. Also, comparing with other existing results on mirror-field entanglement such as |10' we note that in the 
present setup in the absence of a continuum of field modes our results do not indicate the existence of a steady-state 
entanglement at high temperatures. 

While the parameter values chosen here may pertain to a narrow parameter regime corresponding to weak coupling 
and isolation from environment, we have illustrated that there is a significant effect of the idf parameters on the 
mirror-field entanglement. We discuss our results further and conclude in the following section. 


V. DISCUSSION 


The foremost theme in our analysis is to highlight the significance of the internal degrees of freedom of a mirror 
that play the role of the essential intermediary when it comes to studying the interaction between a quantum field 
and the mirror’s mechanical motion. We illustrate how a microscopic model of quantum optomechanics, such as 
the MOF model proposed by Galley, Behunin and Hu |T| is a physically more complete and intuitive description 
for optomechanical interactions, in that not only can it agreeably reproduce the known optomechanical properties 
both in the classical and quantum regimes, it also elucidates new physical aspects which are not accounted for in 
the general description of optomechanical interactions via radiation pressure coupling. Specifically looking at the 
quantum entanglement between the mirror’s mechanical motion and the field, we find that there is a significant and 
even a critical role played by the internal degree of freedom in certain parameter regimes as it can act as a means to 
coherently transfer correlations between the field and the mechanical degree of freedom. 

The MOF model allows us to go beyond the usual disjoint treatment of mirror-field interactions wherein one 
imposes boundary conditions on the field and treats the mechanical effects of the field arising from the radiation 
pressure force separately to attain a self-consistent depiction where we see both the radiation pressure (section IIAI 
and the boundary conditions (section |II B[ ) emerge from a physically motivated microscopic interaction. The new key 
aspects that arise from this self-consistent treatment of the mirror-field interaction can be summarized as follows. 


• Coherent transfer of excitations - We show that the conventional boundary condition approach arises as the 
limiting case of the MOF model where the quantum fluctuations of the idf are lost to the bath. As illustrated in 
Fig^] isolating the idf from the environment can provide an additional channel for coherent transfer of mirror- 
field correlations. Since the zd/-field dynamics is at a much faster time scale as compared to the center of mass 
motion, for an undamped idf we observe much faster time scales for the entanglement dynamics determined by 
the effective id/-field coupling (oof) rather than those from the conventional radiation pressure coupling (~ 13). 

• Fully dynamical description - For relativistically moving mirrors, as in the case of dynamical Casimir effect 
|33] . applying boundary conditions is an inadequate description of the dynamics since in the time scales over 
which internal degrees and the field reach a steady state thereby leading to an effective boundary condition, 
the mirror center of mass moves appreciably enough to affect their interaction. In cases where the timescales 
of the mechanical motion and the field-internal degree of freedom interaction are close to each other, including 
the internal degree of freedom becomes relevant as the only means to capture the coupled dynamical interplay 
of the three subsystems. 


• Field frequency shift - In the MOF description, we observe an additional shift to the field frequency from its 
second order interaction with the idf as seen from (27) and (291, a feature that is not accounted for in the 
boundary condition treatment. Such a diamagnetic term contribution can be significant in the strong coupling 
regimes, even leading to change in the radiation pressure force from attractive to repulsive as has been studied 
in [55] , 


• Radiation pressure shot noise - We also observe that for an undamped idf the radiation pressure force is 
determined not only by the quantum fluctuations of the field but also those of the idf as suggested by ( |40| . While 
in the steady state limit the strength of these fluctuations is largely determined by the boundary conditions, in 
the early time limit the idf being an independent quantum degree of freedom its quantum fluctuations would 
influence the radiation pressure shot noise as well. 


As a result of including an extra quantum degree of freedom one would naturally expect there to be a difference in the 
quantum correlations of the field and the mirror’s mechanical motion. We show that in the parameter regimes where 







16 


the idf is isolated from the environment and for strong coupling, the role of the idf is more pronounced. This can be seen 
from the effective mirror-field interaction strength that is determined strongly by the id/-field resonance condition, 
with the effective bilinear CoM-field coupling going to a maximum at resonance as in (33). The dependence of the 
optomechanical interaction on the «d/-field resonance leading to an enhanced mirror-field entanglement is something 
that can not be captured in the boundary condition treatment of optomechanical interactions. This is seen in the 
plot for MF entanglement as a function of detuning A in Fig(5] 

Also, the time scales for all dynamics, including that of the mirror-field entanglement is largely determined by the 
effective id/-field coupling ao f and the id/-field detuning A as seen from (44 ). One can make the same observation 
from Fig|4]that for an idf isolated from the environment, the entanglement dynamics are at a much faster timescale 
( otoF ) as compared to the boundary condition approach (13). 

Analogous to the case of optomechanical entanglement in a cavity setup, we find that if the field is red-detuned 
with respect to the cavity frequency such that one facilitates the two-mode squeezing interaction between the field 
and the mechanical mode, the entanglement is sustained for longer times as can be seen from Fig. [5] 

From studying the mirror-field entanglement as a function of the various parameters of the model pertaining to 
the idf and otherwise, we find that the presence of the idf can influence the entanglement dynamics to a significant 
extent. Not only does a microscopic model like the MOF model reproduce the known optomechanical properties 
and provide a more self-consistent approach of studying optomechanical interactions, more importantly it leads to 
qualitatively different physics, specifically in the quantum regime. We conclude that the internal degree of freedom 
being the quintessential mediator of quantum correlations between the mirror center of mass and the field, the MOF 
model gives a physically more complete treatment of the mirror-field entanglement. 
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Appendix A: Logarithmic Negativity 


After t = 0 the interaction is turned on and the three subsystems ( mdf idf and the field) begin to interact with each 
other as the reduced density matrices for each of individual becomes a mixed state. The linearity of the interaction 
terms guarantees that the quantum state of the three harmonic oscillators that starts Gaussian remains Gaussian. 
Thus the dynamics of quantum entanglement can be studied by examining the behavior of the quantity s m and 
the logarithmic negativity Emf [38| : 


E = det 


Vmf 


ft 


ifr A ,, 
-M 


Here A4 is the symplectic matrix 1 < 


Emf = max{0, — log 2 2c_} . 

(— i)u y , Vmf is the partial transpose of the covariance matrix 


Vmf = 


f Vm Vmf\ 

V V MF Vp ) 


as defined in @. (c + , c_) is the symplectic spectrum of Vmf PT + (ih/2)M, given by 


c± = 


Z zt \/ Z~ — 4 det Vm f 


1/2 


(Al) 

(A2) 

(A3) 

(A4) 


with 


Z = det Vujif + det Vpj? — 2 det Vmf- (A5) 

For the quantum oscillators in Gaussian state, Emf > 0, E < 0, and c_ < h/ 2, if and only if the quantum state of 
the two subsystems is entangled |39J . Emf is an entanglement monotone m whose value can indicate the degree of 
entanglement. 
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Appendix B: Comparison of the MOF model with boundary conditions 


Let us start with the standard optomechanical treatment where we treat the mirror field interaction via the radiation 
pressure. In 1+1 D, the Hamiltonian of a scalar field interacting (that corresponds to the vector potential of the optical 
field) interacting with a point-like mirror is given by 


Hbc =™(P 2 + Z 2 ) 


dx (^ + J c2 ( dx ®) 2 ~ ^~ d ^ d ^ Z zPMZS(x - Z) 


(Bl) 


Again, assuming that a single field mode is being driven at a frequency ui, we have 

Hbc = ^ (P 2 + Z 2 ) + hcoa^a - h ZzpmAq ) (, ae + aV^ 0 )! 


(B2) 


The pre-factor /?mf = j^ZzpmA oe 1 11 ^ 0 in the interaction term is the standard optomechanical coupling as in |Jj with 


Zz pm = yjjfo as the zero point motion length for the center of mass motion and Aq as the dimensionless field 

amplitude. Moving to a rotating frame with respect to the free field Hamiltonian (Hp = h/jja) a) leads us to the 
following equations of motion for the dimensionless field and mirror variables 


dZ 

d t 

d P 

df 

d$ 

d t 

dn 

d t 


= 13 P 


= — 13Z + \/ 2Ile/+yf p d* — \/ 2 In i/dv/ /. n IT 
= V^Im/?MF-Z' 

= VzR,&I3mfZ 


(B3) 

(B4) 

(B5) 

(B6) 


It can be seen that in the weak coupling limit (A 2 <C 2mujeoc), redefining the field amplitude in terms of the di¬ 
mensionless amplitude A 0 as <h 0 = ij 2 Je 0 L ^°’ the effective mirror-field coupling coefficient in (331 reduces to 


olmf ~ —Pmf/V^- Solving for the entanglement in the two cases, including a damping coefficient for the inter¬ 
nal degree of freedom and setting the idf detuning to be large we see a perfect overlap of the log negativity from either 
approach in Fig|4] 
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